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We investigate the separability properties of quantum two-party Gaussian states in the framework 
of the operator formalism for the density operator. Such states arise as natural generalizations of 
the entangled state originally introduced by Einstein, Podolsky, and Rosen. We present explicit 
forms of separable and nonseparable Gaussian states. 
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I. INTRODUCTION 

In their reasoning concerning the alleged incomplete- 
ness of quantum mechanics |I|, Einstein, Podolsky, and 
Rosen (EPR) used this wave function for a system com- 
posed of two particles: 



*(a;i,X2) 



dpe 



{2ni/h){xi - X2-\- xo)p 



(1) 



It is a singular function of the distance xi — X2 and could 
be visualized as an infinitely sharp Gaussian wave func- 
tion of the entangled two-party system. Bell inequalities 
of some kind are violated for this wave function, as can 
be demonstrated by using its Wigner representation B. 
Recent applications of entangled two-mode squeezed 
states of light for quantum teleportation js) and other 
quantum information purposes [l| have generated a lot 
of interest in the separability properties of general mixed 
Gaussian states in quantum optics [p|. In one approach, 
the separability properties of continuous- variable systems 
in states described by Gaussian Wigner functions have 
been investigated with the aid of Heisenberg uncertainty 
relations 0]. Another approach made use of the crite- 
rion of positivity under partial transposition M. Both 
approaches use the basic definition, namely that a gen- 
eral quantum density operator of a two-party system is 
separable if it is a convex sum of product states [k[ : 



Y.p>^p':'®pf' 



with 2^ Pk — 1 and pk > , 

k 



(2) 



where ph. and p^ are statistical operators of the two 
subsystems in question. 

The authors of H and W arrived at essentially the 
same conclusions, while using very different techniques. 
The objective of this Brief Report is to show how equiv- 
alent results are derived by employing the powerful alge- 
braic methods of quantum optics. We use a direct op- 
erator method to study the separability of an arbitrary 
two-party Gaussian operator G{a, a\ h, b^) of unit trace, 
referring, for instance, to two modes of the radiation field, 
parameterized by their ladder operators a and b. An ex- 
plicit algebraic form of the Gaussian operator enables 



us to decide whether G is a density operator and, if so, 
whether it is separable, that is: whether G is a positive 
operator of the form (|2|). Our method works for arbi- 
trary Gaussian operators , but for the sake of clarity 
and also in view of the importance of the EPR wave func- 
tion (|l|), here we shall carry out the explicit calculations 
only for a specific class of Gaussian operators that form a 
natural generalization of the original EPR state (1^) . We 
wish, however, to stress that the algebraic approach is 
quite general, and that the method provides an explicit 
construction of the Werner decomposition (||) for sepa- 
rable two-party Gaussian states. The algebraic approach 
provides a natural link between the partial-transposition 
criterion of Peres ||l^ and P-representable Gaussian op- 
erators. 



II. GAUSSIAN OPERATORS. BASICS 

Following Wigner flit] , we associate a real phase space 
function W{a, (3) with any such operator, and in partic- 
ular with the Gaussian operators G of interest, in accor- 
dance with fn2l 



G = 2^ d^a d'pW{a,P) 



X : e 



-2{at - a*){a - a)„-2(fet - /3*)(6 - /?) 



(3) 



Here, the integrations are over the two phase spaces of 
the oscillators, parameterized by the complex variables a 
and P, and the normally ordered exponential operators 
: exp(-2(at - a*){a - a)) : , : exp(-2(&t - /3*)(6 - /3)) : 
are parity operators ||l3||, 
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(4) 



displaced in phase space by a and (3. We take for granted 
that the traces of aG and bG vanish; otherwise a unitary 
linear transformation would enforce this condition. 

A Gaussian operator G, then, is one whose Wigner 
function VF is a Gaussian function of its complex vari- 
ables a and f3: 



l^(a,/3) = ^i^e-H^^ 



(5) 



where v^ = [a*,a,P*,P] is a complex 4-componeiit row 
and W > is a positive 4 x 4-niatrix. The positivity of W 
and the prefactor in (0) ensure the correct normahzation 
of G to unit trace, 



Tr G = fd^a f(f(3 W{a, /3) = 1 , 



(6) 



We shall also find it useful to work with the Weyl- 
Wigner characteristic function 

G{a, 13) - Tr {e""^ ~ ^'^Ge^''^ " ^*^} = e^i^^Vv ^ (^^ 

a Gaussian function as well, related to the Wigner func- 
tion W{a, P) by two- fold complex Fourier transforma- 
tion. Accordingly, we have 



V = EW"^E>0 and W = EV"^E, 



(8) 



where E = diag[l, —1, 1, —1] is a diagonal 4x4 matrix. 
Given VF(a,/3) and C(a,/3), Eqs. (|) and (@) do not 
specify W and V uniquely. The symmetry of our stan- 
dard form of V, 
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(9) 



with real ni, n2 and complex mi, m2, Wg, rric, exploits 
this arbitrariness conveniently. 



III. GAUSSIAN OPERATORS. EXPLICIT FORMS 

As stated in the Introduction, we shall illustrate the 
algebraic method by the example of a generalized version 
of the EPR wave function (^ . This generalized Gaussian 
operator is specified by a V matrix of the restricted form 



V = 







n+h 





m* 





m* 

m n + 



m 









n+h 



(10) 



correspondiiig to ni — n2 — n, mi — 7712 = rus = 0, and 
TOc = ?Ti in (y). The constraint 



n + ^ > \m.\ 



(11) 



ensures V > 0. The significance of n and tti, 

n = Tr{ataG} =Tr{6^6G}, m = -Tr{a6G}, (12) 
is revealed upon expanding (uh in powers of v. 



In view of (||) , the matrix W of the Wigner function (|^) 
is at hand, and then (0) gives us the operator in normally 
ordered form. With the identity mm 



-(^a'a 



(i^cr 



(13) 



valid for all complex C (the ( = 2 case is met in (0)), 
we so arrive at one explicit form of the corresponding 
Gaussian operator G, namely 



G = 



(n + l)2-|m|2 
n(n -|- 1) — \m\^ 






■ fotfe 



{n + iy - \m\ 
More compactly, this appears as 

G ^ o Gi G2 o , 



' (n + l)^-|m 



■^ab 



(14) 



(15) 



where the basic Gaussian operators Gi , G2 have the form 



Gi = (l-.gi)<?r" and G2 = (1-52)52 



fotfe 



with 



51 =92 



n{n + 1) — |?np 
{n + 1)2 - |m|2 



and the sandwiching operator S is 



V'(n + l)2-|m|2 



(n + l)2-|m| 



Tab 



(16) 



(17) 



(18) 



Quite generally, the basic Gaussians of (|l^) have a fi- 
nite trace if —1 < 51,52 < 1- More specifically, the con- 
straint (11) implies here that — (4n-|- 3)~^ < 51 = 52 < 1 
and thus ensures the finite value of the trace. Positiv- 
ity of Gi and G2 , and therefore also of G itself, requires 
more restrictively that < 51 , 52 < 1 in ( [iq ) , irrespective 
of the particular form that S might have. For V of the 
specific form (|l^), this says that G > is equivalent to 



n{n + 1) > \m\ or n > y |m|2 + 



As a compact statement about V, this appears as 



V 



iE>0 



(19) 



(20) 



with the diagonal matrix E of (g). 

Note that, if (pi] ) is obeyed but (|l^ ) is not, then we 
have a Gaussian operator that does not represent a den- 
sity operator although its Wigner function is positive and 
properly normalized because the matrices V and W are 
positive. These matters are illustrated in Fig. m 

The explicit construction of the Gaussian operator was 
here performed for a matrix V of the specific form (M). 
In the most general situation of (0), we have more pa- 
rameters, but G is always of the generic form (Iq), that 
is: a product G1G2 of two thermal Gaussian operators, 
sandwiched by a 5*^, S pair. 
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FIG. 1. Concerning the parameters of the Gaussian opera- 
tor associated with the V matrix of (|lO|). Only n, |r7i| values 
above the dashed line are allowed by constraint ( |ll| ) . Accord- 
ing to ([19|), values on or above the solid line specify positive 
Gaussians of the form (p3). For values on the solid line, the 
Gaussian operator is a projector. Separable Gaussians belong 
to values on or above the dash-dotted line; see (p2). 
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INOPA) =c"'('^^^^^"'')|0,0) = — ^(tanhr)"|n,ri), 



T!=0 



(25) 



where A = tanh r relates the squeezing parameter r to A 
and thus to parameter n of the Gaussian operator. 

We said above that (|lj) is a natural generaUzation of 
the EPR state (|^). The stage is novif set for justifying 
this remark. To this end, we first note that, for real A, 
the position wave function of |*) of (p3) is given by 



*(a::i,a;2) = ^= exp 

\/7r 



(1 + X'^)ixl + x\) - 4Aa;ia::2 



2(1 -A2) 

— - /dpe^*/^)^^^ -X2)p 
nh J 



x,/lH?e-^[(P/")' + i(^i+^^)'l, (26) 
1 + A 



and then observe that 



(27) 



obtains in the limit A ^ 1. Indeed, this is the EPR state 
of dil) with xo = 0. 



IV. GAUSSIAN OPERATORS. PURE STATES 



V. SEPARABILITY OF GAUSSIAN STATES 



In this section we investigate Gaussian operators that 
are projectors and thus represent pure states. This case 
occurs when the equal sign holds in (|l9|), so that gi = 
(72 — 0, and 

(l-5i)5f'^-0'^''^ = <5. as 51-0, (21) 



for example, states that 
G1G2 



0"''^ + '''=^'5.t,,o'5,t,,o = |O,O)(O,0| (22) 



in this limit. The Hilbert space vector \n,m), here for 
n = TO = 0, denotes the state with n quanta of a-type 
and m of 6-type. 

As a consequence, Eq. ( [l4|) turns into 



G= |*)(*| = (1 - |Ap)e^"^''V:0)(0:0|e 



X*ab 



(23) 



where A — —m/{n +1). Such an operator projects onto 

oo 



(24) 



n=0 



For real A, we recognize an example of the well known 
two-mode squeezed state that can be generated by Non- 
degenerate Optical Parametric Amplification, 



A positive Gaussian operator G is said to be P-repre- 
sentable if it can be written in the following form: 

G= fd^a fd'^(3P{a,P) 

X : e-^*^^ - °*)('^ - °)e-(^^ " ^*)('' " ^) : , (28) 

with a non-negative phase-space function P{a, P) that 
must not be more singular than a Dirac i5 function. 
The ordered exponentials : exp(— (a^ — a*){a — a)) : and 
: cxp(— (6^^ — /3*){b — /?)) : are projectors onto the coher- 
ent states labeled by a and /3, respectively. 

For a P-representable Gaussian operator, we have 



P(«,^) = ^l|Ze-HPv 



(29) 



with P related to the 4x4 matrix V of the characteristic 
function (R) and the matrix W of the Wigner function 
©by 



E(V-il) ^E= (W-i 



I 



(30) 



where I is the 4x4 unit matrix. So, a given Gaussian 
operator is P-representable if 



V-il>0 



(31) 



If the left-hand side is truly positive, we have a four- 
dimensional Gaussian in (^^ , else it is the product of a 
two-dimensional Gaussian and a two-dimensional S func- 
tion, or the product of two two-dimensional S functions. 
For the Gaussian operator (|lj), the existence of the 
P-rcpresentation is guaranteed if 



n > \m\ . 



(32) 



In Fig. n these n, \m\ values are on or above the dash- 
dotted line. Since the ordered exponentials in (ES) 
project to coherent states, such P-representable Gaus- 
sians are convex sums of product states. They are thus 
of the separable kind, as defined in (0), where the for- 
mal summation over k is now the two-fold phase-space 
integration of (P8|). 

As Peres observed flfl], the partial transpose p^" of 
a separable statistical operator (0) is another statistical 

(k) 

operator because the pa 's are turned into other positive 
operators and the p^ s are not affected to begin with. 
In other words, p"^" > is a necessary property of a sep- 
arable p. As surmised by Peres and demonstrated by the 
Horodecki family p5[ , it is in fact sufficient for two-party 
systems composed of two spin-i objects ("qubits") or of 
one qubit and one spin-1 object. 

Concerning the systems of interest here, of two har- 
monic oscillators, the Peres criterion p^" > does not 
imply that p is separable. But, as Simon noted [Q, in 
the particular case that p is a positive Gaussian operator, 
Peres' criterion is sufficient to ensure that p is separable. 
Indeed, for a positive Gaussian G the Peres criterion is 
equivalent to requiring that G is P-representable. 

To see this, let us be more specific and agree on using 
the Fock representation to define the partial transpose. 
Then, G -^ G^" amounts to VK(a,/3) -> W{a*,(3) in §), 
that is: 



W^T„,WT„, with T„ 
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(33) 



in (H) and, as @ implies, V -^ ETaE VETaE in (g). 
For the V matrix of ([lO|), this results in 



G' 



1 



n + m 



i{at-^6t)(a_^*6) 



n + 1 + |r7i| \n + 1 + |?7i| 

1 ( n— |m| \ 2 



n-\-\ — \m\ \n-\-\ — \m\ 



i(at+/xfet)(a + ;,*6) 



(34) 



where [i = ml\m\. Therefore, we have G"^" > only if 
n > |rn|, so that G is not separable for n < \m\. In view 
of (p2h, then, the partial transpose is positive whenever 
the Gaussian in question is P-representable. And, as 



already remarked after (|3^), G is separable if it is P- 
representable. Together these observations say this: 

A positive Gaussian operator is separable 

if it is P-representable, and only then. (35) 

This statement is more generally true than our argument 
suggests, because the limitations that originate in the 
special form of V of (IG) can be lifted 0] . 

Note, in particular, that the projectors of Sec. IV are 
non-separable for \ra\ = ^n(n + 1) > which includes 
the EPR limit of n — > oo. The Gaussian projector (ESI) is 
separable only in the other limit of m = 0, n = 0, when 
it projects onto the two-oscillator ground state |0,0). 
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